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We point out mistakes made in the one-loop calculation of some diagrams for the process ny* — > 
7 in the preprint arXiv:0807.0296 , and present correct results. Especially, we have difficulty to 
understand their argument that a highly off-shell gluon generates a light-cone (infrared) singularity. 
It is shown by means of the Ward identity that the gauge-dependent light-cone singularity found in 
larXiv: 0807.0296 does not exist. It is then shown that a hard kernel derived in the kr factorization 
of exclusive processes is gauge invariant and free of the light-cone singularity. 

o : 

' In a recent preprint larXiv:0807.0296l [l| , the authors studied the pion transition form factor in the fcy factorization 
theorem at one-loop level, calculating the diagrams in full QCD for the form factor (see their Fig. 4) and the diagrams 
for the fcy-dependent pion wave function (see their Figs. 2 and 3) in the covariant gauge. They focused on a special 
piece of contribution, which depends on the gauge parameter and contains a light-cone singularity. It was found that 
<^ the full QCD diagrams do not generate this contribution, but those for the pion wave function do, if the involved 
_ partons are taken to be off-shell. Therefore, the hard kernel, defined as the difference of the above two sets of 
diagrams, is gauge dependent, and contains the light-cone singularity. They were then led to the conclusion that the 
kr factorization theorem for exclusive processes violates gauge invariance, and that the perturbative QCD (PQCD) 
approach 0, i, i, [| to exclusive B meson decays, based on the kx factorization theorem, is also gauge dependent. 



I This contradicts the conclusions drawn in our previous work Q. In this comment we point out the mistakes made 
Q-i^ in their calculation, and demonstrate that the gauge-dependent light-cone singularity discussed by the authors of [l| 
^ ■ does not exist. The all-order proof for the gauge invariance and infrared finiteness of a hard kernel derived in the kx 
i ' factorization ^ is indeed correct. 

Consider the process 7r(P)7* —> y{p) in Fig. [U with P (p) being the pion (outgoing on-shell photon) momentum 
' along the plus (minus) direction. Define the momenta of the valence quark and anti-quark in the pion as [l| 

^: fc5' = (fc+,o,fcii) , ^ik+,o,-h±) , (1) 

respectively, with k^ = xqP'^ and fc^ = (1 — xo)P^ , xq being the momentum fraction. The leading-order (LO) hard 
kernel is given, in terms of the above momenta, by 



00 

o ; ^(0) ^ I ^ I ^2) 

g . 2ktp- + kl^ xoQ2 + kl^ ' 



\ with the momentum transfer — 2P~^p~. At next-to- leading order (NLO), a loop momentum q carried by the 
' ^ ' additional gluon may flow through the hard kernel, for example, in Fig. [21 In this case the factorized fcr-dependent 
' wave function is convoluted with 
H ■ 

■ ■ ■ ■ = -— . (3) 

2(fc+ - g+)p- + Ifcii - gl|2 

We find that the mistakes made in the calculation of Fig. [5] and Fig. [3] 1] arise from an improper application of the 
contour integration. Take the simple loop integral associated with Figs. [2l^d), [Sj^b) and[3l^e) as an example: 

/d^Q 1 
(27r)4 {q^ - XI + i£)[q^ - (1 - a)Xl + le] ' ^"^^ 

where the gauge parameter a comes from the gluon propagator in the covariant gauge [l| 



g2-(l-a)Ai+: 



(5) 
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FIG. 2: NLO diagrams for the pion wave function with the loop momentum flowing through the hard kernel. 



and the gluon mass Xl is introduced to isolate the infrared (IR) pole. Applying the Feynman parametrization, it is 
easy to obtain 

Aryn/ A? 

/=^ln^+UVpole, (6) 

where InA^ denotes the IR singularity [l|, the explicit expression for the ultraviolet (UV) pole is irrelevant here, and 
II is the renormalization scale. However, if employing the contour integration naively, Eq. (j4]) vanishes: 

r^-lf,. 2 f I (7\ 

^'"^^ J (27r)4 (2q+q- - ql - Xl + ^e)[2q+q- - ql - {I - a)Xl + ^e] " ' ^ 

The poles in the above integrand are located in the same half plane, either upper or lower, no matter whether one 
integrates over the q^ or q^ first. Therefore, one can close the contour in the other half plane not containing the 
poles. Obviously, Eq. ([7]) contradicts Eq. ([6]). How to understand this puzzle is the key to clarify the conflict between 
Ul and @. 

The above simple example illustrates that a naive application of the contour integration in the light-cone coordinates 
may lead to a false result. This is exactly the reason the authors of [l| drew the wrong conclusion from Fig. [21 When 
one closes the contour in Eq. ^ by including the semicircle at infinity, it has been implicitly assumed that the 
contribution from the semicircle is negligible. However, this assumption does not hold as g"*" — > (3]: the product 
q^q^ in the denominator does not provide a suppression on the semicircle of a large radius |g~| as g"*" — > 0. In other 
words, the poles q~ = {q\ + X\)/{2q'^) and q~ = [q\ + (1 — a)A|]/(2g+), also moving to infinity as q'^ 0, cannot 
be avoided by the contour in Eq. ([7]). A safe way to proceed is to close the contour of q^ with a semicircle of a 
large but finite radius R in the half plane without the poles. For example, the semicircle in the upper half plane is 
considered for (7+ > in Fig.dl^a). The integration over q~ along the real axis is thus equated to the integration over 
the semicircle, and Eq. ([7]) becomes 



/ — 16iag'^ lim 
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de dq+ + i dd dq^ 







— (8) 



{2q+Re'^ -q\-Xl+ ie)[2q+ Re'" -q\-{l- a)A| + ie] 



where the first and second terms in the above square brackets correspond to Fig. IHa) and Fig. |D^b), respectively. 
Performing the integration over g"*" and 9, and then taking the R 00 limit, we reproduce Eq. ([5]). With the above 
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FIG. 3: NLO diagrams for the pion wave function without the loop momentum flowing through the hard kernel. 
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FIG. 4: Contours for Eqs. (0 and dH) in the regions (a) > and (b) < . 



prescription, the potential contribution from the large semicircles at g+ — > is taken into account, and a consistent 
result is obtained in the contour integration. 

Below we recalculate Fig. [5] (the NLO diagrams for the pion wave function with the loop momentum flowing through 
the hard kernel) following the aforementioned method, and demonstrate that Figs. [2l^a), [2fb), and [2jc) arc in fact 
free of the gauge-dependent light-cone singularity alnXj^. We also recalculate Fig. [3] (the NLO diagrams for the pion 
wave function without the loop momentum flowing through the hard kernel), and point out the erroneous observation 
derived from Figs.[3ljc) andlS^f) in Since the authors of [l| have evaluated the NLO diagrams for the form factor 
using the Feynman parametrization, these results are valid. 

We start from Fig.^^b), whose gauge-dependent part is written as [l| 

J. I iR- 2 /■ d.'^l 2 (fc+ - g+)q- - • + , , +^^r2/r ft w fn\ 

= mag^ J (2^)4 [(fc,-g)2+,,](,2+,,)2 ^(fc^ - (fcf - {k^ (fcr^ - q^)) . (9) 

The wave function (j)a\2b is convoluted with the LO hard kernel: 



Ca 2b ' 



>H^°'> = f dx ( d^k^ / ,2 0a|2b , (10) 



with the variable x = /P~^. Integrating the two (5- functions over x and k±, the LO hard kernel in Eq. ([3|) appears. 
In the light-cone region the scaling law for the components of q is defined by {q~^ ,q~ , q±) ^ {S'^,1,S) with 6 being a 
small parameter according to which the first term {kf — q^)q^ in Eq. ^ gives the leading contribution. The 
light-cone singularity is regularized by introducing the gluon mass as shown in Eq. ([S]). Applying the contour 



(a) (b) (c) 

FIG. 5: Contours for Eqs. ^ and (I12|l in the regions (a) q+ > k+, (b) < g+ < k+, and (c) q+ < 0. 



integration over q for < 9+ < naively, we reproduce the result in 

<^r^|2. ® H("^ = ^^iL + finite terms . (11) 

P+TT XqQ'' + kf^ 

As explained above, we have to examine the potential contribution from the large semicircles. For the range 
> k'^ , all the poles are located in the lower half plane, so we close the contour of q~ with a large semicircle in the 
upper half plane shown in Fig. El^a). The contribution from this semicircle remains negligible in the limit g"*" — > fcj*", 
because of the suppression factor {k^ — q'^) in the numerator of Eq. ([9]). For g+ < 0, all the poles are located in the 
upper half plane, so we close the contour with a large semicircle in the lower plane shown in Fig. [S{c). In this case, 
the semicircle may contribute as — > 0, and we equate the integral over q~ to the contribution from the semicircle 
of a radius R. For < (7+ < , the single pole q~ = {\ki±_ — — *^)/[2(<7''' — ^1")] is located in the upper half 
plane and the double pole q^ = {q"^ — ie)/(2q+) in the lower half plane. We close the contour of q~ through the upper 
half plane with a semicircle shown in Fig. [5]Jb), which picks up the residue from the single pole, giving Eq. pT|) . The 
contribution from the semicircle of a radius R needs to be subtracted, since it does not vanish at infinity as q~^ — > 0. 

Therefore, Fig. EJb), as convoluted with the LO hard kernel, contains the additional contribution from the two 
semicircles, one for < q^ < k^ and another for q^ < 0, 

i dd dq+ + i dd dq+ 

J TV Jo J — TT J —00 

2{k+ - q+){Re''>f 

(2g+i?e'» -ql- \l)[2q+ Re'<^ - ql - {1 - a)\l][2{q+ - k+)Re'0 - \kii_ - giP] 



^ (fc+ - q+)QyP+ + - giP ■ 
Working out the integration over q^ , 9 and q±, and then taking the i? — > cx) limit, we obtain 

(I}'j2b® H<-°> = - ^2 + finite terms , (13) 

which cancels the IR pole in Eq. pTjl exactly, implying that Fig. [DJb) is free of the gauge-dependent light-cone 
singularity. Note that the limit i? — ^ 00 must be taken after the q~^ integration. Then the contribution from the 
semicircles at (7+ = can be evaluated correctly. The same conclusion applies to Fig.[2I^c). 

We have explained the absence of the light-cone singularity in Fig. [2lb) through the contour integration over q~ 
by including the semicircle contribution. Each of the residue contribution in Eq. (jlip and the semicircle contribution 
in Eq. has the light-cone singularity, but they cancel each other exactly. It is thus understood why the authors 
in [l[ got a fake light-cone singularity from Fig. [l)Jb): they have missed the contribution from the large semicircles. 
Their wrong result in Eq. (jll[) has an origin similar to that of Eq. ([7]). Performing the integration over q~^ first, 
we obtain the same result as in the case starting with the q^ integration. In the contour integration over g"*", the 
semicircle contribution must be handled appropriately as q^ 0. For the range q^ > 0, we equate the integral 
over q~^ to the contribution from the semicircle of a radius R in the upper half plane. For q~ < 0, we close the 
contour of g"*" through the lower half plane with a semicircle of a radius R, which picks up the residue from the pole 
g+ — {k^Q'^ + P^\ki± —q±\'^ — ie) /Q^. It is found that neither the residue contribution nor the semicircle contribution 
contains the light-cone singularity, because the integral in Eq. Q may be singular in the region (g"*", q^ ,q±) ^ (5'^, 1, (5), 
but not as g~ — > 0. 



,/ I tt(q) 16iag^ f d^q 
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If applying the contour integration naively to Fig. ^c) , 



(0) _ 



P+ 



2(k^ - q^)q - ki± ■ q± 



9i 



(2^)4 [(fci - 9)2 + is]{q^ + ze)2(xoQ2 + fc2^) 



(14) 



the mistake made in [T| is more obvious. The term {k^ — q^)q~ gives a light-cone singularity from the region 
{q^ , q~ ,q±) ^ (5^,1,(5) in the naive contour integration, which is the same as in Eq. Ijlip . Since the loop momen- 
tum does not flow through the hard kernel, the term (k^ — q'^)q^ also generates a UV divergence from the region 
{q~^,q~,q±) ^ (1,A^,A) with the scale A ^ oo in the naive contour integration. However, according to the covari- 
ance argument [l|, the corresponding loop integral, proportional to q~ , should vanish like k^ = 0. To overcome 
this apparent contradiction, the authors in [l[ dropped the gluon mass Al, adopted the dimensional regularization, 
misinterpreted the UV divergence as another light-cone singularity, and made the UV divergence and the light-cone 
singularity cancel each other. The momentum configuration q (1, A^, A) should give a UV divergence, since it must 
be regularized with the number of dimensions n < 4 in the dimensional regularization, and cannot be regularized by 
the gluon mass A^. Hence, we have difhculty to understand their argument [l| that a gluon with infinite invariant 
mass 5^ A^ ^ cx) produces the light-cone (infrared) singularity. The fact is that Fig. [3jc) has neither UV divergence 
nor light-cone singularity, after carefully including the contribution from the semicircles. 

Our result is natural from the viewpoint of the Ward identity. Contracting the loop momentum q to the vertex on 
the internal quark line (see the numerator of the gauge-dependent term in Eq. ([5])), 



3- /:i 



A-, 



5- ^1+ /l i>- Jzx 



3- 



(15) 



(p - k^f " (p - fci + qf {p - fci)2 {p-ki+ qf ' 

the first (second) term leads to the LO hard kernel associated with Fig. [3l^c) (Fig. [2jb)). The diagram Fig. 4c in [l[ 
for the form factor is then factorized into convolutions of the LO hard kernel with Figs. (H^b) and[31Jc). If Fig. 4c in 
[H and Fig. M.c) do not contain the light-cone singularity as observed in [l|. Fig. [^Jb) should not either. Another 
support to our result comes from the following simple observation. In the light-cone region the q^ and dependence 
in the LO hard kernel is negligible. Then Figs.[2lb) andlS^c) involve exactly the same loop integral, namely, the same 
behavior in the light-cone region. If Fig. does not produce the light-cone singularity as found in 1], Fig. Hfb) 
should not either. 

We then turn to Fig. , which appears as a consequence of factorizing the box diagram Fig. 4a of [1] for the 
form factor. Again, Fig. [2fa) should not contain the gauge-dependent light-cone singularity, because Fig. 4a of [l| 
does not. The explicit expression for the leading contribution from Fig. [IJa) in the light-cone region is given by 



{kt-q+)(k+ + q+){q-y 



(2^)4 [(fci - 9)2 + i£][{k2 + qY + ie\{q^ + ief 

Xd{k+ - {k+ - q+))S\k^ - (h± - q±)) , (16) 

where the numerator ((7~)^ arises from q'^q'^ in the gluon propagator for p = u = — . Integrating over q^ first, we 
have two poles, q^ = {\kij_ - (LlP - ie)/[2{q+ - fcj*")] for < q+ < k^ and q^ = {\kij_ - (f_Lp - ie)/[2{q+ + fc^)] for 
— fc^ < g"*" < 0. Closing the contour in the upper half plane for the former, and in the lower half plane for the latter 
naively, the result in [1] is reproduced. 



4Q;a, 



InAi 



P+TT XoQ^ + kjj 



finite terms 



(17) 



which is also a fake singularity. 

Similarly, the contribution from the two semicircles of a radius R, which were added to form the closed contours 
mentioned above, is written as 



6^|2a®H^°^ = ^ lim 



P+ R^oo J {2tt 



d^qi_ 
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ie\d, 



-2{k+ - q+)Re^o - Ifcij. - qi 



2{k+ + q+)Re^^ - |fci_L - + is 



1 



{2q+Re^ -q\-\l + ie){2q+Re'0 - ql - {I - a)\l + ie) 
1 



(fc+-g+)g2/p+ + |fcii-9l|2 



(18) 
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The rest of the procedure is straightforward, leading to 

The sum of Eqs. ([TT]) and i.e., Fig. [2I^a), is free of the gauge-dependent hght-cone singularity, contrary to the 

observation in [1]. In summary, all the diagrams in Figs. [2] and [3l except those with the gluons attaching only the 
Wilson lines (Figs. Hfd), [3]^b), and[3]^e)), do not generate the gauge-dependent IR singularity. 

We then comment on the diagrams Figs. [2l^d), [3Kb), and[3Ke). These diagrams do not appear in the factorization of 
collinear gluons [1, [l^ from the pion transition form factor. It has been shown that their sum is IR finite [1]. This 
must be the case, since all the IR divergences in the form factor diagrams, which arise from the collinear region with 
the scaling law {q~^,q~,q±) ~ (l,<5^,i5), have been absorbed into the other diagrams in Figs. [2] and [3] [6]. However, 
in view of a gauge invariant definition for the pion wave function, Figs. [IKd), [S^b), and[3Ke) should be included. A 
resolution is to invoke a soft subtraction factor in the denominator 

Six -Xk )^ i '^^y^ ^-i(l-x)P+v--ik, -gi (0|g(y)W^yW^-^n;y,oWo(n)7+75g(0)|g(fci)g(fc2)) 
J (2^)2 {mv{nVWy{u)In-,yMU,W,{n)W,{um ' 

which removes Figs. HJd), [3Kb), and [3Ke) in a gauge invariant way. In the above definition, y = (0,y~,yj_) is the 
coordinate of the anti-quark field q , n with ^ the direction of the Wilson line, and \q{ki)q{k2)) the leading Fock 
state of the pion. The factor Wy (n) denotes the Wilson line operator 



Wy (n) — P exp 



igs I dXn ■ A{y + An) 







(21) 



The two Wilson lines Wy{n) and Wo{n) must be connected by a link In.y.o at infinity (Tol. [Tl|. 

The subtraction factor in the denominator generates the soft diagrams similar to Figs. [2| and [3l but with the fermion 
lines being replaced by the Wilson lines in the direction u. Besides Figs. [21(d), [S^b), and[3Ke), additional soft diagrams, 
such as the vertex corrections with the gluons attaching the Wilson lines in the directions n and u, are introduced at 
the same time. To ensure that the soft subtraction does not change the collinear structure of the numerator, u should 
not lie on the light cone, namely, ^ 0. Other than this requirement, the direction of u is completely arbitrary. 
Hence, the subtraction of Figs, [l^d), [Sib), and[3Ke) will result in an arbitrary UV pole for the pion wave function. 
One can then take advantage of this arbitrariness, adjusting u so that the UV pole of the above vertex corrections 
cancels the UV pole of the self-energy corrections to the Wilson lines along u. Below we demonstrate this cancellation 
explicitly. Because the subtraction factor is gauge invariant, we evaluate the soft diagrams in the Feynman gauge, 
i.e., consider only the g^'^ tensor for the gluon propagator. Then the former is given by jl2i] 



a^Cp u-n , \/ (u ■ n)^ — v?n? -\- u ■ n /I , Ann? . 

\ L . . T+lnTT^ , (22) 



and the latter by 



asC^fl.^^P^. (23) 
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It is easy to find that the above two expressions cancel exactly for u^n'^ < and u-n/ ^\u?'n^\ « 0.85. It is also trivial 
to show that the remaining soft diagrams with the loop momentum flowing through the hard kernel cancel, when the 
above conditions for u and n are satisfied. Eventually, we just need to calculate the seven effective diagrams in Fig. 2 
of 0, i.e., Figs.[2Ka)-[2Kc), and Figs.[3Ka),[3Kc),[3Kd) and[3Kf) in this comment, in the Feynman gauge at one-loop level. 

After clarifying the conflict, the all-order proof of the leading-power kT factorization theorem for the pion transition 
form factor in Sec. HI of [6] follows straightforwardly. Below we mention the points of the proof briefly. Since the 
light-cone singularity does not exist, the substitution for the metric tensor of the gluon propagator in the covariant 
gauge, Eq. (47) of [6], works for extracting all IR (collinear) divergences. We then employ the derivative in Eq. (50) 
of [6] to collect the gauge-dependent terms arising from the off-shell partons. The derivative of the soft subtraction 
factor vanishes, since it is gauge invariant. The derivative of the pion wave function is then related to the derivative 
of the numerator in Eq. (|^ . which, after applying the Ward identity, gives Eq. (53) in Combining Eqs. (52) and 
(53) in the fc-r factorization theorem for the pion transition form factor is proved by induction. 

In conclusion, the only necessary revision for [6| is to replace the definition for the fcj^-dependent pion wave function 
by Eq. ([^0)1 with an appropriate vector u. Except this modification, all the results in Q are valid. The gauge- 
dependent light-cone singularity found in [l[ does not exist; their result is attributed to a careless application of 
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the contour integration in the hght-cone coordinates. The correct prescription is to keep the contribution from the 
semicircles of a finite radius first, work out the loop integration, and then move the semicircles to infinity. We have 
confirmed that the all-order proof of the fcy factorization theorem for the pion transition form factor holds: the gauge 
dependence, arising from the off-shell external partons, cancels between the full QCD diagrams for the form factor 
and the effective diagrams for the pion wave function. Therefore, the kx factorization produces a gauge invariant and 
IR finite hard kernel, and provides a solid basis for the PQCD approach to exclusive B meson decays. 
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